Spectral-fluctuations test of the quark-model baryon spectrum 
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We study the low-lying baryon spectrum (up to 2.2 GeV) provided by experiments and different 
quark models using statistical tools which allow to postulate the existence of missing levels in 
spectra. We confirm that the experimental spectrum is compatible with Random Matrix Theory 
-the paradigmatic model of quantum chaos-, and we find that the quark models are more similar 
to a Poisson distribution, which is not compatible with what should be expected in a correlated 
spectrum. From our analysis it stems that the spectral fluctuation properties of quark-model spectra 
are incompatible with experimental data. This result can be used to enlighten the problem of missing 
resonances. 
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Since the discovery of the first excited state of the 
nucleon, the A(1232) baryons have played a cen- 
tral role in the study of the strong interaction. After 
this discovery, baryons proliferated and their mass spec- 
trum attracted a lot of attention. In the 60's and 70's, 
nonrelativistic quark models were worked out 0] and 
the quark model framework was established, contribut- 
ing strongly to the development of Quantum Chromo- 
dynamics (QCD). Later on, quark models evolved into 
relativistic versions 0, 0, S| • 

It is well known that the number of baryons predicted 
by quark models 0, 0, 0, 0] is substantially larger than 
what is observed in meson scattering and production ex- 
periments Q . This fact raises the problem of missing res- 
onances, which has opened the door to a huge experimen- 
tal effort in the last years to observe and identify these 
missing states p}. These experiments have to achieve 
high precision due to the important background (which 
can veil resonances) and the overlap of baryons, as well as 
the need to survey different meson production channels 
and observables. The procedure to assess the existence of 
these elusive baryons consists on analyses of partial waves 
[H and polarization observables of the reactions com- 
paring experimental data from different sources to what 
is obtained after including or removing the hypothetical 
resonance. If data are better reproduced, the existence 
of a resonance is possible but sometimes debatable (the 
pentaquark search is a very good example of how difficult 
these studies are [HI). The Particle Data Group (PDG) 
rates the possible existence of the resonances based on 
the quantity and quality of experimental data. Only after 
several independent experiments and analyses, a baryon 
is awarded a well-established status, rated with three or 
four stars. 

In the last years, Lattice QCD (LQC D) flll| and effec- 
tive QCD-inspired models (EQCDiM) [12| have started 
to provide insight on the non-perturbative energy regime 
of QCD where the low-lying baryons live, but they are 
still far away from providing a complete analysis. Hence, 
we have to resort to tractable models, such as quark mod- 



els, to study the baryon spectrum. For this reason, quark 
models are expected to retain in forthcoming years the 
importance they had in the study of the static properties 
and internal structure of baryons. 

It is well established that baryons are aggregates of 
partons. Consequently, the mass spectrum of low- lying 
baryons is an energy spectrum of a standard many-body 
quantum system, like an atomic nucleus, and it consists 
on all the possible bound and excited states which stem 
from an interacting many-body quantum system. Since 
Wigner discovered that the statistical properties of com- 
plex nuclear spectra are well described by Random Ma- 
trix Theory (RMT) (l3j |. statistical methods have become 
a powerful tool to study the energy spectra of quantum 
systems. The most striking result in this field is that the 
statistical properties of the energy-level fluctuations are 
universal and determine if a system is chaotic or inte- 
grable. Integrable systems display a non-correlated se- 
quence of levels, which follows the Poisson distribution 
14j , whereas chaotic systems are characterized by a cor- 
relation structure described by RMT [15]. This kind of 
analysis has been already applied to the hadron mass 
spectrum by Pascalutsa [l6| obtaining a chaotic-like be- 
havior. 

In this Letter we apply the spectral statistical tech- 
niques to the problem of missing resonances. Recent 
works by Bohigas and Pato [T3] prove that if we ran- 
domly remove some energy levels from a correlated spec- 
trum, it partially loses correlations and becomes closer 
to a Poisson distribution regardless its actual correlation 
structure. Following this result, we can obtain hints on 
the existence of missing resonances comparing the spec- 
tral correlations of the baryon spectra supplied by ex- 
periments and quark models. We have analyzed the ex- 
perimental spectrum (Breit- Wigner masses) provided by 
the PDG (set EXP in what follows) Q and the spectra 
given by the relativistic quark models by Capstick and 
Isgur (set CI) Q and by Loring et al. [f| (sets LI and 
L2 which correspond, respectively, to models A and B 
in H). We have considered all the resonance states in 
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these spectra up to 2.2 GeV. In principle, if the exper- 
imental spectrum is not complete due to the absence of 
some baryons which are not observed but are predicted 
by quark models, and if we assume that missing reso- 
nances are randomly distributed [HI, the experimental 
spectrum should be less correlated than the theoretical 
ones. 

Prior to any statistical analysis we have to accomplish 
two preliminary tasks. First of all, it is necessary to iden- 
tify the different symmetries involved in the spectra. If 
a sequence of levels involves more than one symmetry, 
its spectral statistics are deflected towards a Poisson dis- 
tribution (see HI, IH for generic reviews and for a 
recent work where the effects of both mixing symmetries 
and missing levels in the same sequence are surveyed). 
Hence, it is necessary to extract from the full spectrum 
sequences of levels involving the same symmetries (quan- 
tum numbers) to proceed with the spectral analysis. 

In the spectra considered in this Letter, we can iden- 
tify the following symmetries associated to the baryons: 
spin, isospin, parity, and strangeness. Strangeness can be 
dropped due to the assumption of flavor SU(3) invari- 
ance. Therefore, for every statistical analysis, we split 
all the spectra in sequences where all the levels present 
the same values of spin, isospin, and parity. For reasons 
stated below we only account for sequences with three or 
more levels. 

The second preliminary task is the unfolding proce- 
dure. In any energy-level spectrum, we can split the level 
density, p(E), into a smooth part, ~p{E), and a fluctuat- 
ing part, p(E), with p{E) = p(E) + p(E). The unfold- 
ing procedure allows to extract the fluctuating part from 
the level density, removing the smooth component of the 
spectrum. There are several ways to unfold a spectrum 
and we choose the simplest one. First, we compute the 
distance between two consecutive levels, Si = Ei + \ — Ei, 
and then we rescale Si using its average value Si = Si/ (S) 
[2l| . The resulting quantities are called Nearest Neighbor 
Spacings (NNS). This procedure undergoes some prob- 
lems, specially in the long-range correlation analysis [22| . 
but it is suitable for the kind of analysis of our concern 
in this Letter. 

From the NNS we obtain one of the most relevant 
quantities in spectral statistical analysis: the Nearest 
Neighbor Spacing Distribution (NNSD). The NNSD fol- 
lows the Poisson distribution P(s) = exp(— s) if the spec- 
trum is integrable (non-correlated) [14|, but it follows 

the Wigner surmise P(s) — ^ exp ( —iI f~) ; which stems 
from RMT, for a chaotic (correlated) spectrum [l5|- For 
our purpose here, it is enough to consider that the less 
(more) correlated the sequence of levels is, the closer to 
the Poisson (Wigner) distribution the NNSD is. 

In order to obtain a significative result we have calcu- 
lated the NNSD, P(s), for each one of the four different 
spectra we survey: EXP, CI, LI, and L2. We account for 



all the sequences {si}xi where X stands for the quantum 
numbers which identify each sequence [23[ . Set EXP has 
70 energy levels distributed in 15 sequences; set CI, 145 
levels and 19 sequences; set LI, 142 levels and 21 se- 
quences; and set L2, 104 levels and 19 sequences. 
We also have evaluated the function 
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which is related to the accumulated NNSD and allows a 
better study of the tail of the distribution. 

Results are depicted in Fig. [T] The experimental spec- 
trum is closer to the Wigner surmise than to the Poisson 
distributionfin good agreement with what Pascalutsa 
obtains in [16]). On the other hand, the three theoretical 
spectra are quite close to the Poisson distribution, so that 
they are less correlated than the experimental spectrum. 
These results are opposite to what is expected from an 
experimental spectrum spoiled by missing resonances. If 
the experimental spectrum is not complete due to miss- 
ing states, it has to be closer to the Poisson distribution 
than the theoretical ones. 

However, there is one subtle objection to our analysis 
regarding the unfolding procedure which has to be con- 
sidered in detail. As we have pointed before, the unfold- 
ing has to be performed independently for each {si}x 
sequence. This entails that, at least in some cases, we 
have worked with very short sequences of levels. In such 
cases, the unfolding procedure could give rise to mislead- 
ing results if it makes the spacings Si spuriously closer 
to Si ~ 1, shifting the NNSD far away from the Pois- 
son distribution and bringing it to the Wigner surmise 
(roughly centered at s = 1). As an extreme example, let 
us see what happens to a sequence of two levels after the 
unfolding. In this case, we only have S\ — E2 — E\, and, 
therefore, s% = Si/ (S) = 1. This is the reason why we 
have included only sequences with three or more levels in 
our analysis. Moreover, since the experimental sequences 
are shorter than the theoretical ones, they can be more 
affected by this problem, spoiling the direct comparison 
we have made above. Therefore, to reach a final conclu- 
sion it is mandatory to perform an improved analysis. 

To elucidate whether the unfolding procedure can give 
rise to misleading results, we have accomplished a more 
sophisticated analysis of the NNSD. Since the unfolding 
could yield different effects on experimental and theoret- 
ical spectra, we have avoided a direct comparison among 
them. Instead, we have built four Poisson-like and four 
RMT-like reference spectra, each of them optimized to 
study each one of the four sets of data (EXP, CI, LI, and 
L2). We have proceeded in the following way for every 
set. For each sequence {si}x of a given spectrum we have 
built two equivalent reference sequences {s^}x with the 
same length, one RMT-like and another Poisson-like. In 
doing so, we start from one long RMT-like spectrum and 
one long Poisson-like spectrum. We divide each spectrum 
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TABLE I: Probability to obtain, under the null hypothesis, a 
value of the Kolmogorov-Smirnov test statistic as extreme as 
that observed. 



(a) Experimental values from PDG (set EXP). 




(b) Model by Capstick and Isgur (set CI). 




(c) Model by Loring et al. (set LI). 




(d) Model by Loring et al. (set L2). 

FIG. 1: NNSD for the experimental spectrum provided by 
the PDG data, the model by Capstick and Isgur, and two 
parametrizations of the model by Loring et al. The histogram 
represents the spacings; the solid line, the Wigner surmise; 
and the dashed line, the Poisson distribution. The inset shows 
the function F(x) in logarithmic scale. 



in as many subsequences as sets with the same quantum 
numbers the spectrum under consideration has, each one 
with its appropiate length. Next, we unfold each ref- 
erence sequence independently and calculate the NNSD 
for each set. Finally, we compare each spectrum ob- 
tained directly from the experimental data (EXP) and 
the quark models (CI, LI, and L2) to the two ad hoc 
built-up RMT-like and Poisson-like reference spectra. In 



Spectrum 


EXP 


CI 


LI 


L2 


Poisson 


0.25 


0.44 


0.37 


0.10 


RMT 


0.56 


4.5 ■ 10 -6 


1.3 ■ 10 -4 


6.5 • 10" 3 



this way, each reference spectra is distorded by the un- 
folding in the same way as sets EXP, CI, LI, and L2 are. 
Consequently, we can verify whether their correlations 
are true or spurious due to the unfolding. 

In order to obtain significative results, we have per- 
formed the Kolmogorov-Smirnov distribution test 24| in 
each case. As null hypothesis we have chosen that the 
studied NNSD coincides with the ad hoc built-up refer- 
ence Poisson-like or RMT-like distributions, against the 
hypothesis that both distributions are different. In Ta- 
ble U we summarize our results. We observe that all the 
spectra emerging from the three theoretical models are 
incompatible with RMT, whereas the experimental one 
seems to be closer to the Wigner surmise than to the Pois- 
son distribution. Set L2 provides a result that seems to 
be incompatible with both Poisson and RMT behaviors. 
From panel (d) in Fig. [T] it stems that the probability 
to observe spacings closer to zero are higher than for a 
Poisson distribution. Hence, the deviation from a Wigner 
surmise is even larger than from a Poisson distribution. 

These results, together with those shown in Fig. [TJ 
point out that the experimental spectrum is more corre- 
lated than the three theoretical ones. Following the work 
of Bohigas and Pato [ItJ , we can say that this is incom- 
patible with the usual statement that there are missing 
resonances in the experimental spectrum which are in- 
cluded in the theoretical models. Moreover, the Hamil- 
tonians used in such models do not describe the statistical 
properties of the experimental spectrum; they rather cor- 
respond to an integrable system whilst the experimental 
spectrum is close to a chaotic system. Hence, quark mod- 
els, as they are presently built, may not be suitable to 
reproduce the low-lying baryon spectrum, and, therefore, 
to predict the existence of missing resonances. 

On the other hand, it is important to notice that quark 
models assume triality (three-quark states), while QCD 
allows non-three-quark states such as pentaquarks and 
hybrid baryons (excited glue). Therefore, according to 
QCD degrees of freedom, the spectrum which is obtained 
within quark models cannot be complete. Consequently, 
the NNSD of any quark model (whatever the interacting 
Hamiltonian is) may be deflected to a Poisson distribu- 
tion due to its inherent uncomplctcncss. The importance 
and quantification of such effect depends on the amount 
of missing non-three-quark states present in the actual 
spectrum and remains as an open question. However, we 
cannot blame the non- Wigner character obtained for the 
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quark models studied in this Letter on this effect: the 
loss of correlation cannot be due to missing non-three- 
quark states because theoretical models predict more lev- 
els than what is observed. 

Regardless we obtain that present quark models are 
not able to predict missing resonances in the experimen- 
tal spectrum, our results are compatible with the exis- 
tence of some missing states. In fact, the shape of the 
NNSD and the Kolmogorov-Smirnov test for a RMT-likc 
spectrum with a 20% of random missing levels arc alike 
to those of set EXP. 

The analysis presented in this Letter should be ex- 
tended to baryon spectra provided by LQCD and 
EQCDiM as soon as complete calculations become avail- 
able. In this way we can use the universal properties of 
fluctuations to examine the results given by these models 
or to use the spectra provided by LQCD and EQCDiM 
to test the universality of fluctuations. 
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